Confinement effects on glass forming liquids probed by DMA 
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Many molecular glass forming liquids show a shift of the glass transition T g to lower temperatures 
when the liquid is confined into mesoporous host matrices. Two contrary explanations for this effect 
are given in literature: First, confinement induced acceleration of the dynamics of the molecules 
leads to an effective downshift of T g increasing with decreasing pore size. Secondly, due to thermal 
mismatch between the liquid and the surrounding host matrix, negative pressure develops inside 
the pores with decreasing temperature, which also shifts T g to lower temperatures. Here we present 
novel dynamic mechanical analysis measurements of the glass forming liquid salol in Vycor and 
Gelsil with pore sizes of d — 2.6, 5.0 and 7.5 nm. The dynamic complex elastic susceptibility data 
can be consistently described with the assumption of two relaxation processes inside the pores: A 
surface induced slowed down relaxation due to interaction with rough pore interfaces and a second 
relaxation within the core of the pores. This core relaxation time is reduced with decreasing pore 
size d, leading to a downshift of T g oc 1/d in perfect agreement with recent DSC measurements. 
Thermal expansion measurements of empty and salol filled mesoporous samples revealed that the 
contribution of negative pressure to the downshift of T g is small (< 30%) and the main effect is due 
to the suppression of dynamically correlated regions of size £ when the pore size d approaches £. 

PACS numbers: 64.70.Pf, 61.20.Lc, 62.25.+g 



I. INTRODUCTION 



When approaching a glass transition some physical prop- 
erties like viscosity or relaxation times change up to 
14 orders of magnitude 1,2 . An explanation for the ob- 
served slowing down of the dynamics is the formation 
of collectively dynamically rearranging clusters^ or re- 
gions, with growing size £ and increasing relaxation 
times as T g is approached 5 . The idea of an increas- 
ing dynamic correlation length £ when approaching a 
glass transition is strongly supported by recent computer 
simulations^—, although not strictly proven, since com- 
puter simulations cannot treat the time range of the a- 
process. Very recently a breakthrough was achieved in 
this field. G. Biroli, et al£ found first-time evidence 
that the Mode Coupling Theory predicts a growing dy- 
namic length scale approaching the glass transition of 
a supercooled liquid. The authors obtained a rather 
modest growth of the dynamical length scale £ with 
decreasing temperature, which is in very good agree- 
ment with computer simulations^ and experimental re- 
sults. Indeed many experimental setups like heat capac- 
ity spectroscop y 10 i 11 , multidimensional NMRi 2 -^ 4 , mul- 
tipoint dynamical susceptibilities 15 , etc. have been used 
to monitor a possible growing length scale accompany- 
ing the glass transition. All these results agree in the 
fact, that the obtained dynamically correlated regions - 
although material dependent - are of the order of 1-4 nm 
and display - if at all - a weak temperature dependence. 

An alternative experimental approach to get a refer- 
ence to a possibly existing cooperation length £ which 
increases when T — »• T g is by spatial limitation of a 
glass forming liquid. Spatially confining geometries as 
ultrathin films, mesoporous silica or zeolithes have al- 
ready been used to study phase transitions of water 16 , 



o 20 



9 

an 



hydrocarbons , noble gase o i , liquid crystals^ or 
alkenes 21 . But this concept also illuminated the old and 
still open question on the very nature of the glass tran- 
sition and its dynamics 22 : In a pioneering work Jack- 
son and McKenna 23 studied the glass transition of or- 
ganic liquids in controlled pore glasses (CPG) for vari- 
ous pore sizes d. They found a reduction of the glass 
transition temperature T g for liquids in confinement as 
compared to the bulk material. The downshift of T, 
was larger for smaller pore sizes, i.e. AT g oc 1/d, 
effect similar - but not as strong as - the suppres- 
sion known for the melting temperature T m in confine- 
ment. During the following two decades this effect was 
studied via calorimetry 24 - 25 , dielectric spectroscopy^ 6 -, 
neutron scattering 2 ^, light scattering 2 ^ and molecular 
dynamics 2 ^. It was shown that in many cases confine- 
ment below a characteristic length impedes 24 the transi- 
tion, implying that molecules within a region of the size 
(approaching T g typically some n m 30 i 31 ) have to co- 
operate and rearrange in order to create the glassy state. 
Hindering this cooperation first leads to a downshift of 
T g if d ~ t; g and finally to a suppression of the transition 
if d < 

However, although this shift of T g with decreasing 
confinement size was found in many studies, there are 
complications which blur this simple picture: E.g. in 
some systems a competition appears between slow- 
ing down of molecular motions due to pinning of the 
molecules at the pore surface and acceleration of the dy- 
namics due to decreasing size of the confinement. An- 
other effect occurs due to the difference in thermal ex- 
pansion coefficients of the porous host matrix and the 
glass forming liquid. This may create negative pressure 
upon the confined liquid when the glass transition is ap- 
proached. Some authors attribute the whole observed 
downshift of T g to this negative pressure effect 28 . We 



will address these points in very detail below. For excel- 
lent reviews about these topics the reader is referred to 
Refs. [Handlll 

Very recently the confinement effect on the glass- 
forming liquid salol was studied via dynamic mechanic 
analysis measurements (DMA) 34 in Vycor with d = 7 
nm pore size. It turned out that the dynamic elastic 
response is very sensitive to the glass transition of li- 
quids confined to mesoporous samples. Based on the re- 
sults of computer simulation s 29 ! 35 we could disentangle 
acceleration effects due to confinement and slowing down 
of molecular motion due to interaction of the molecules 
with the rough pore surface. We could even predict the 
pore size dependence of the dynamic elastic response (see 
Fig. 4 of Ref. [34J). In order to test these predictions and 
to study the glass transition of salol for different pore 
sizes, further measurements have been performed. Here 
we present novel experimental results of the tempera- 
ture and frequency dependence of the complex dynamic 
elastic susceptibility of salol confined in mesoporous ma- 
trices of d = 7.5 nm, 5.0 nm and 2.6 nm. In addition, 
thermal expansion measurements have been performed, 
which now allows us to take a new look at the often dis- 
cussed negative pressure effect on glass forming liquids in 
confinement and to seperate this effect from an intrinsic 
size effect. 

The present paper is organized as follows: Section II 
yields insight into sample preparation and some technical 
details of DMA analysis. Section III displays a compila- 
tion of the experimental data and results of modeling 
and interpretation of the present data. It also contains 
a calculation of the effect of adsorption swelling and the 
separation of the actual downshift of T g in salol into the 
negative pressure effect and the confinement effect. Sec- 
tion IV concludes the paper. 



II. EXPERIMENTAL 

A. Sample preparation 

Porous Vycor samples are made by Corning Inc., NY, and 
sold under the brand name "Vycor 7930". Via phase sep- 
aration and leaching a three dimensional random network 
of pores in nearly pure silica is fabricated 3 ^. Pores are 
uniformly distributed in length, direction and density 37 . 
The mean ratio of average pore diameter d and pore 
length I, is d/l « 0.23. Gelsil is a mesoporous xerogel 
consisting of pure silica with a very narrow pore radius 
distribution. Gelsil rods were made by 4F International 
Co., Gainesville, FL. Results on pore sizes derive from 
BJH analysis of the individual N2-desorption isotherms 38 
and are summarized in Tab. HI 

All samples were cut and sanded in order to gain 
parallel surface plains. The typical size of a sample was 
(2x2x8) mm 3 for parallel plate and about (2x1x7) mm 3 
for three point bending DM A- measurements. Cleaning 
was done in a 30% hydrogen peroxide solution at 90° C 
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Tab. I: N2 adsorption characteristics of porous silica samples. 





Gelsil2.6 


Gelsil5 


Vycor 


av. pore diameter (nm) 


2.6 


5.0 


7.5 


surface area (m 2 /g) 


586 


509 


72 


pore volume (cm 3 /g) 


0.376 


0.678 


0.214 


porosity <p 


0.51 


0.66 


0.30 




Fig. 1: Sketch of (a) parallel plate and (b) three point bending 
geometry. 



for 24 h, drying at 120° C in a high- vacuum chamber 
at 10 -6 bar, also for about 24 h. The guest glass form- 
ing material was salol (phenyl salicylate, C13H10O3), a 
low molecular weight liquid, whose melting temperature 
is T m = 316 K. Salol is a standard, so called fragile 39 , 
glass former (m=73) known^ 5 - to form a glass either at 
extreme cooling rates of 500 K/min or in pores smaller 
than 11.8 nm. Filling was done at 317 K via capillar- 
ity wetting. By comparing the weight of clean and filled 
samples the filling fractions / were determined (see Table 

in). 



B. Dynamic mechanical analysis (DMA) 

In this method a static and a dynamic force F stat + Fd yn ■ 
e iut (0.001-16 N at 0.01-100 Hz) are applied on a sample 
using a quartz or steel rod (see Fig. Q}. The response of 
the sample is measured via the displacement of the rod. 
Absolute height h, height amplitude Ah and phase lag 5 
are read via electromagnetic inductive coupling (LVDT) 
with a resolution of 10 nm and 0.01° respectively. These 
data allow direct access to real and imaginary part of 
the complex elastic susceptibility at low frequency and 
as a function of temperature and applied force. In ad- 
dition, the thermal expansion of a sample can be deter- 
mined in the so called TMA-mode, where no external 
force is applied. Two devices are used: a DMA 7 and a 
Diamond DMA, both from Per kin Elmer Inc. Two mea- 
suring geometries are applied: parallel plate compression 
(PP) and three-point bending (3PB), see Fig.Q] 

Parallel plate geometry reveals purely the complex 
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Young's modulus Y* = Y' + iY" ', where Y 7 and Y" are 
the storage and the loss modulus, respectively. The three 
point bending geometry delivers Young's modulus plus 
a small (geometry dependent) contribution of a shear 
elastic constant. More details on measurement geome- 
try may be found in Refs. H3 and The absolute ac- 
curacy of resulting real and imaginary parts Y' and Y" 
is rather poor, mainly because of contact losses between 
the quartz rod and the sample. A discussion of these 
systematic errors may be found in Ref. B In contrast 
the relative accuracy is excellent and the DMA- method is 
estimated to be about 100 times more sensitive to detect 
glass transitions or other subtle phase transitions than 
DSC-measurements 4 ^. 



III. RESULTS AND DISCUSSION 
A. Dynamic elastic response 

Diamond DMA-measurements (in parallel plate and 
three-point-bending geometry) of Vycor and Gelsil sam- 
ples filled with salol are shown in Figs. [2]- [5] The loss mo- 
dulus Y" (Fig.[2p) of salol in 7.5 nm pores clearly shows 
a "two peak structure", i.e. a peak with HWHM about 
20 K, and a shoulder or second peak at about 15 K higher 
temperature (also see Fig. |5p). This is also reflected by 
the real part Y ; , which displays a "two step like shape" 
with temperature (Figs. and[5£i). Both peaks in Y" 
shift to higher temperatures with increasing frequency as 
expected for a glass transition. In smaller pores of Gelsil 
5.0, peak and shoulder merge into one asymmetric peak 
of width ~ 30 K (see Fig. and [5^), also shifted with 
higher frequency to higher temperatures. In 2.6 nm pores 
the loss peak shows a rather symmetric form broadened 
up to about 60 K, see Figs. Hp and [5f. 

While in large pores of 7.5 nm diameter vitrifica- 
tion of salol seems to happen decoupled (two peaks in 
Y") in regions near the pore surface and the pore center, 
things change in smaller pores. With decreasing pore 
diameter, Y" approaches a symmetric form and simul- 
taneously Y' changes from a "double step shape" into 
a "single step shape", indicating only one type of relax- 
ation process. Similar broadening effects as for the loss 
peaks of our DMA-measurements were observed in pores 
of decreasing size also by calorimetric 23 and dielectric 
measurements 4 ^ 4 -. This broadening as well as a shift of 
the glass transition to lower temperatures was calculated 
by Sappelt and Jackie using kinetic Ising and lattice gas 
models^ 5 ., and shown to originate from confinement in- 
duced suppression of cooperative motion of molecules. 

Pure Vycor and Gelsil, meaning exposed to air and 
therefore mostly filled with nitrogen, do not show any 
of these features. Y' decreases about 2% between 300 
K and 180 K. Y" is constant within the corresponding 
temperature range. 

Any standard relaxation model like Debye, Kohlrausch, 
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Fig. 2: Real (a) and imaginary parts (b) of the complex 
Young's modulus of Vycor 7.5 nm filled with salol (filling 
fraction / « 0.79) measured in three point bending geometry. 
The curves are offset from the 1 Hz data for sake of clarity. 



Cole-Cole or Cole-Davidson fails to describe our dynamic 
elastic susceptibility data if only one type of relaxation 
process is assumed. One would have to use extreme 
stretching parameters to fit Y', which then leads to im- 
proper temperature shifts of the peaks in Y" with respect 
to the experimental data and misfitting signal heights. 
The most efficient model to describe our data turned 
out to be a modification of the empirical Vogel-Fulcher- 
Tammann-law 



t(T) = r • exp[ 



E 



T-Tn 



(1) 



where To is a preexponential factor, E-ku is an activation 
energy, and T is the Vogel-Fulcher (VF) temperature. 
Following computer simulation s 29 ! 35 we take into account 
a shift of VF-temperatures along the pore radius r . In 
a recent paper Zorn et al. 46 suggest the empirical ansatz 



T (r) = Too + 



k 



R - r + r v 



(2) 



with the bulk VF-temperature Too, and the pore radius 
R = d/2. The so called penetration radius r p is the radius 
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Fig. 3: Real (a) and imaginary (b) parts of the complex 
Young's modulus of Gelsil 5.0 nm filled with salol (filling 
fraction / « 0.75) measured in parallel plate geometry (Dia- 
mond DMA). 1Hz signal are original data, other signals are 
offset for sake of clarity. 
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Fig. 4: Real (a) and imaginary (b) parts of the complex 
Young's modulus of Gelsil 2.6 nm filled with salol (filling 
fraction / « 0.65) measured in three point bending geometry 
(Diamond DMA). 1 Hz signal are original data, other signals 
are offset for sake of clarity. 



beyond which it is very unlikely to find a particle in the 
fluid stated 9 . The combination of Eqns. JT} and (0 leads 
to a radial distribution of relaxation times r inside the 



pore: 



r(r, T) = r • exp \ 



E 



T - (Too + 



(3) 



R-r-\ 



Eqn. (0 describes the exponential increase of relaxation 
time when a rough pore wall is approached, and a grow- 
ing influence of the pore wall with decreasing tempera- 
ture, a behaviour which was also found by recent com- 
puter simulations^^^. A Cole-Cole-Plot of log(Y') 
vs. logiY") of our data clearly calls for a Cole-Davidson 
model of the complex dynamic elastic susceptibility 



Y*(u) oc 



1 



(1 + icur) 2 



(4) 



with u = 27tz/, v being the measurement frequency, and 
the broadening parameter 7. Using Eqn. J3]), averaging 
over the pore radius R, and separating real and imaginary 
part of Y* = Y' + i • Y" by common procedures leads to 

2 Ay f R cos[y arctan(cor(r,T))] 
Y ~ 1 " — Jo [1 + ^,^]? ^ (5a) 
2 AY f R sinb-arctanjuT^m 
Y ~ — Jo [l + ^Mr.D 2 ]* 1 ' 

As already mentioned above, the two peak structure in 
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Y" of 7.5 nm and 5 nm confined salol (Figs. [2] and [3]) sug- 
gests to split the dynamic elastic response into a core and 
a surface contribution: the molecules in the core (center 
of the pores) behave bulk-like and are dynamically de- 
coupled from the molecules near the pore surface. This 
is modelled by inserting into Eqn. ([5b]) the corresponding 
relaxation times To • exp(E / (T — To)) given by Eqn. (p} if 
r < R c and r(r, T) given by Eqn. J3]) if r > R c (see also 
Fig. 7). The sum of the two contributions perfectly de- 
scribes our Y' and Y" data on salol in 7.5 nm and 5 nm 
pores simultaneously (see Fig.0. 

In 2.6 nm pores we do not expect any molecule to 
behave like the bulk liquid any more, since the pore ra- 
dius is of the same order as the estimated surface shell 
(see Tab. ITT|) , implying that every molecule is influenced 
by the near surface. Therefore we use Eqns. (5) with no 
bulk term which reproduces one single peak and also fits 
our data very well (Fig. [5fc and f and Fig. [6b and f). 

The radius of the "core" of bulk- like interacting 
molecules turned out to be R c = 2.5 nm and 1.35 nm 
in 7.5 nm and 5.0 nm pores, respectively, see Tab. [Til 
This implies that the thickness of the shell of molecules 
being slowed down by wall interaction R — R c = 1.25 nm 
and 1.15 nm for 7.5 nm and 5.0 nm pores respectively. 

Additional loss peaks, attributed to molecules form- 
ing H-bonds to the inner pore surface, have also been 
reported from dielectric measurements of salol in 7.5 nm 
pore o 48 i 49 . The work of Kremer and Stannarius^ 9 - also 
revealed that the typical size of a shell of molecules in- 
teracting with the pore surface is about 2 or 3 molecules. 
Since the size of a salol molecule is estimated as (1.4 x 



0.6 x 0.4) nm 3 in Ref.^or as 0.282 nm 3 in Ref.l5li both 
corresponding to a mean diameter of 0.8 nm, this shell 
size is in the order of 1.6 to 2.4 nm. This is in very good 
agreement with our findings (see Tab. HI]). The core size 
R c decreases with decreasing pore radius (see Tab .Inland 
Fig. O, also in very good agreement with the results of 
Kremer et al. 

Tab. II: Fit parameters used in Eqns. (5) for fits of Fig. [6j 





Vycor 


Gelsil5 


Gelsil2.6 


R (nm) 


3.75 


2.50 


1.28 


r v (nm) 


0.36 


0.25 


0.28 


E(K) 


1750 


1750 


1750 


Too (K) 


158.5 


156.2 


136.0 


r (s) 


10" 11 


10" 11 


10" 11 


7 


0.33 


0.18 


0.15 


k (nm- K) 


18 


11 


25 


R c (nm) 


2.5 


1.35 




shell R - R c (nm) 


1.25 


1.15 


1.28 



The fitted Vogel-Fulcher temperature T 00 is reduced with 
respect to the bulk and with decreasing pore size (see 
Tab. HII). I n order to compare our results with published 
data, we plotted the relaxation time in the pore center 
r(r = 0, T) for various pore sizes d and determined the 
corresponding T g (d) by using the common procedure^ 
for finding the laboratory glass transition temperature, 
i.e. a cut with a horizontal line at r = 100 s (see Fig. [8]). 
As shown in Fig. [TTJ this leads to glass transition tern- 
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Fig. 7: Modelled relaxation time distributions in pores of diameter 7.5 nm to 2.6 nm from Eqn. ((3]) used in Eqns. (5) for fits 
of Fig. m 
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Fig. 8: Relaxation time in pore centers calculated from 
Eqns. ((TJ) and ((3j) with corresponding parameters from Tab.lTTl 
Horizontal line shows r = 100 s. 
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Fig. 9: Height of Vycor sample during the filling process 
against y/i. Inset shows sample height against time. 



peratures decreasing oc 1/d in very good agreement with 
published data of DSC measurements 28 . 



B. Filling process and accompanying effects 

By using a DMA in a static TMA mode one can detect 
small changes in a sample's height with a resolution of 
10 nm. We measured the time dependent swelling of the 
Vycor and Gelsil samples during filling with salol and 
the thermal expansion of empty and filled samples in the 
following way: In parallel plate mode, the quartz rod 
is placed on top of the sample with force F = N, and 
just height and temperature signals are read out. A clean 
piece of Vycor /Gelsil sample is cooled down to 170 K. Af- 
terwards the sample is heated slightly above the melting 
temperature T m = 316 K of salol and kept there isother- 
mally. Crystalline, powder-like salol placed right around 
the sample melts and percolates the Vycor/Gelsil sam- 
ple due to capillarity (Fig. . After filling until satura- 



tion, the sample is cooled down to 170 K again (Fig. [TP]) . 
The time-dependence of the filling process is displayed in 
Fig. [9] for Vycor. While salol is percolating the sample, 
the temperature is held constant and the sample 's height 
is measured. Charts for Gelsil 5.0 nm and Gelsil 2.6 nm 
look very similar. The diagrams in all cases show the typ- 
ical v^-behavior as expected for a single capillary rise ex- 
periment, following Lucas 5 ^ and Washburn 54 . This result 
is in concordance with findings of P. Huber et al*** 5 -, who 
investigated the mass uptake of porous silica samples and 
its time dependence, leading to the Lucas- Washburn y/t 
behaviour of the mass uptake with time. Very recently 
it was shown that the Lucas- Washburn equation (with 
small modifications) works well even at the nanoscale 56 , 
which is in harmony with our results. 

The expansion of a porous sample during adsorption 
of gases or water has already been investigated in the 
1920s 57 . As a liquid/gas intrudes the sample it is sub- 
ject to a negative hydrostatic pressure inside the pores, 
which leads to an expansion of the porous sample dur- 
ing adsorption of gases or water. Mesoporous media 
have enormous inner surfaces up to some 100 m 2 /g (see 
Tab. |TJ). This leads to a considerable stress reduction 
within the whole matrix and a sudden voluminal growth, 
which slows down and stops as all pore space is filled 
(see Fig.0. The change in height due to the adsorption 
swelling can even be calculated quantitatively. The pres- 
sure reduction of the liquid in a capillary is known^ 8 - as 
P c = 2cr/r, with the surface tension a and the capillary 
radius r. With a = 1.73 • 10" 2 N/m from Ref. Hi this 
yields a capillary pressure of 26.6 MPa in 2.6 nm pores. 
This would lead to a hypothetical capillary rise of 1.8 km 
for salol. The linear strain e = Ah/h accompanying the 
filling process can be computed by the equation 5 ^: 



f-Pc 



1 

K 



1 



(6) 



with the filling fraction f, the bulk modulus K of the 
empty host matrix, and the bulk modulus of the mate- 
rial building the solid frame K s (which is nearly pure 
S1O2). The bulk moduli K have been determined by res- 
onant ultrasound spectroscopy (RUS) 60 . Tab. IIHI shows 
parameters used to calculate e = Ah/h. The calculated 
values for the adsorption swelling agree rather well with 
the experimental results. 



C. Negative pressure effect 

The downshift of the glass transition in nm-confining 
pores is often reported to obey a 1/d law, see Refs. [23|, 
H3, HH and H At first this was proposed by Jackson and 
McKenna32, following their former results on the shift 
of the melting transition T m in confinement^ 2 -. But the 
supposed suppression of molecular cooperation when the 
pore diameter approaches an inherent length scale is not 
the only possible reason for a downwards shift of T g in 
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Tab. Ill: Variables of Eqn. ©. 





Vycor 


vjeisno 


vjeisiiz.o 


d (nm) 


7.5 


5.0 


2.6 


porosity <3? 


0.31 


0.66 


0.51 


P c (MPa) 


9.2 


13.8 


26.6 


K (GPa) 


8.1 


3.9 


9.6 


/ 


0.77 


0.62 


0.32 


tcalc 


2.3 • 10~ 4 


6.6 • 10~ 4 


2.2 • 10" 4 


Cexp 


1.0 • 10" 4 


4.1 • 10 -4 


3.5 • 10" 4 



confinement. Zhang et al. proposed the increase of ne- 
gative hydrostatic pressure within the pores due to mis- 
matching thermal expansions of liquid and host matrix 
as the main driving force for the downshift of T g . This 
idea is also discussed by Patkowski et al*2£ and Simon et 
al*££, and was reviewed by Alcoutlabi and McKenna^. 

As Fig. [TUk shows, for large pores and in a cooling 
process starting at RT, at higher temperatures the Vy- 
cor matrix is not affected by its filling. It contracts like 
the empty Vycor matrix with a thermal expansion coef- 
ficient a = Ah/(h ■ AT) = 5.1 • 10" 6 K" 1 . Patkowski 
q\2&, proposed the possible flow and equilibration of the 
confined liquid well above T g , which we also consider to 
be the case here. But as vitrification sets in at about 
230 K, the filled Vycor matrix is subject to a contraction 
which is stronger compared to the empty Vycor sample. 
Strong interaction (H-bondings) between salol and the 
pore surface might be the reason for this. At smaller 
pores of filled Gelsil samples (Fig. [TUb and c) additional 
contraction already starts at higher temperatures. For 
an estimation of the process developing negative pres- 
sure upon the filling liquid, the strain misfit between the 
glass and the host matrix is 

Ae mf = 3(ai - a 2 )AT (7) 

with oti, the thermal expansion coefficients of the host 
matrix (1) and salol (2). Negative pressure then derives 
from AP = Ac/kt, where kt is the bulk compressibility 
of salol. The resulting shift of T g , i.e. 



Results of these estimations are given in Tab. [TVl Pa- 
rameters used for salol are k,t = 5 • 10~ 10 Pa -1 from 
Ref. 0, the thermal expansion coefficient 71 = 3a 1 = 
7.36-10" 4 K" 1 from Ref. 1, and dT g /dP = 0.204 K/MPa 
from Ref. [65|. Our measurements are in accordance with 
enthalpy recovery results of S. L. Simon et al^. Their 
model shows that effective negative pressure develops 2 
to 2.5 K above the reduced glass transition for samples 
with 11.6 nm and 25.5 nm pore sizes. Further, they state 
"...if negative pressure were the cause of the depressed 
T g , the temperature at which isochoric conditions are 
imposed would have to be ~ 20 to 40 K above T g ." For 
comparison we obtain a necessary AT = 10 to 40 K for 
d = 7.5 to 2.6 nm pores, which is in very good agreement 
with Simon et al. 

Tab. IV: Parameters of AT g estimations, AT^ xp = AT™ P + 





Vycor 


Gelsil5.0 


Gelsil2.6 


d (nm) 


7.5 


5.0 


2.6 


a 2 (K" 1 ) 


2.1 • 10" 5 


4.6 • 10" 5 


1.1 • 10" 5 


AT (K) 


6 


10 


15 


Ae mf (%) 


-0.40 


-0.60 


-1.06 


AP (MPa) 
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AT™ P (K) 
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-4.3 


AT conf ( K ) 


-1.4 


-2.9 


-6.4 


ATg Xp (K) 


-3.0 


-5.3 


-10.7 



In our opinion our calculated AT™ P is still overes- 
timated for two reasons: First, using the bulk value a 2 
of the host matrix from Fig. [TU] does not take into ac- 
count internal pore walls being affected by the negative 
pressure inside, relaxing to some extend and so reduc- 
ing pressure. Secondly, thermal expansions of other glass 
forming liquids e.g. toluene have been reported 1.5 times 
smaller in confinement 66 compared to bulk. Moreover, 
thermal expansion of liquid salol dropsi to a quarter of 
its value at the glass transition. So, as the glass transi- 
tion sets in, a\ starts to decrease and a purely pressure 
induced downshift AT g would be even more diminished. 
Apart from this the reason for the size dependence of the 
thermal mismatch effect (see Fig. [TTJ open circles) is not 
clear at all. 



T g (P) = T g (P = 0) 



dTg 

dP 



P=0 



■AP 



(8) 



crucially depends on the choice of AT, the temperature 
range, in which the effective negative pressure upon sa- 
lol develops. This effective temperature range can be 
estimated from our data as follows: As calculated from 
Eqn. J6]) the host porous matrix expands with filling due 
to the negative capillary pressure which acts on the con- 
fined liquid. Since with cooling the liquid salol contracts, 
this stress relaxes and the composite is stress free if the 
filled sample height is the same as for the empty matrix 
which occurs at T = T* (see Fig. So AT « T* - T g . 



IV. CONCLUSIONS 

The glass transition of salol confined to porous host ma- 
trices of Vycor and Gelsil with pore sizes of 7.5, 5.0 
and 2.6 nm has been measured for the first time by 
Dynamic Mechanical Analyzers (DMA 7 and Diamond 
DMA, Per kin Elmer). The dynamic complex elastic sus- 
ceptibility data can well be fitted assuming two types 
of dynamic processes: A "bulk" relaxation in the core 
of the pores and a radially increasing "surface relax- 
ation" of molecules near the pore surface. The calcu- 
lated core relaxation time shows a typical Vogel-Fulcher 
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■ Patkowski et al. 
* present results 

o negative pressure 
effect 



0.1 0.2 0.3 0.4 

-1 -1 
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Fig. 11: Shift of glass transition temperature against 
(pore diameter) -1 . Boxes are T g s from Fig.[HJ triangles show 
literature values from Ref. |28|, open circles display the maxi- 
mum negative pressure contribution (see chapter III. C. ). 



temperature dependence and decreases with decreasing 
pore size d. This confinement induced acceleration of 
dynamics leads to a shift of the glass transition tempera- 
ture T g oc 1/d, which is in perfect agreement with recent 
DSC results 2 ^. Measurements of the sample height with 
filling (adsorption swelling) and thermal expansion are 
used to calculate the effect of "negative pressure" due to 
thermal mismatch between the porous host matrix and 
the glass forming liquid. Such negative pressure could at 
least partly explain a shift of T g in confined glass form- 
ing liquid o 28 i 32 i 63 . Our data show that for salol this effect 
of thermal mismatch could describe at most 30% of the 
observed downshift of T g , which is in harmony with en- 
thalpy recovery experiments^ 3 -. 

In our opinion the main cause for the shift of T g is 
a hindering of cooperativity due to confinement. This 
is also supported by an estimation of this effect using 
the results of Hunt et al. 67 . They calculated the finite 
size effect of the glass transition from percolation and 
effective medium models, which yields 



Inserting t = 100 s, v p h = 1/to, and our fit parame- 
ters from Tab. [HI and assuming that the typical distance 
between molecules ro is about the diameter of a salol 
molecule 51 (d ~ 0.8 nm), we obtain AT^ unt as 3.2, 4.8 
and 9.1 K for 7.5, 5.0 and 2.6 nm pores, respectively. 
These calculated values agree surprisingly well with the 
measured confinement induced downshifts of T g (d) (see 
Fig. EG and AT^ xp in Tab. HV|> . 

Moreover Eqn. J9j predicts 67 that the size depen- 
dence of AT g increases with increasing fragility 68 



In (10) (T g - T ) 2 ' 



since m oc E. Indeed, this correlation between AT g (d) ~ 
m was verified experimentally for many systems, i.e. for 
glycerol 61 (m = 53) AT g {d = 2.5 nm) « -4 K, 
benzyl-alcohol 23 (m = 65) AT g (d = 2.5 nm) « -9 
K, salol 28 (m = 73)AT 5 (d = 2.5 nm) « -11 K, o- 
terphenyi^m = 81) AT g (d = 2.5 nm) w -25 K. 

We think that these considerations, i.e. the down- 
shift of T g calculated via percolation theory, as well as 
the clear correlation between the magnitude of induced 
T g shift and the fragility of a glass forming liquid, both 
confirm our other findings (see Fig. [TTj) that the main 
effect of the confinement is to suppress cooperative mo- 
tion. Negative pressure effects although always present 
contribute only little. 
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